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Dark matter (comprising a quarter of the Universe) is usually assumed to be due to 
one and only one weakly interacting particle which is neutral and absolutely stable. 
We consider the possibility that there are several coexisting dark-matter particles, 
and explore in some detail the generic case where there are two. We discuss how the 
second dark-matter particle may relax the severe constraints on the parameter space 
of the Minimal Supersymmetric Standard Model (MSSM), as well as other verifiable 
predictions in both direct and indirect search experiments. 
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Dark matter (DM) is at the heart of any study regarding the interface between parti- 
cle physics, astrophysics, and cosmology. Its relic abundance has now been measured with 
precision. Combining the results of the WMAP Collaboration and the Sloan Digital Sky 
Survey, QcDAih"^ = 0.110 ± 0.013 (2(t) [1], where ^cdm is the DM energy density normal- 
ized by the critical density of the Universe and h = 0.71 ± 0.05 (2cr) is the scaled Hubble 
parameter. Many dark-matter candidates have been suggested in various models beyond 
the Standard Model (SM) of particle physics, but a nearly universal implicit assumption is 
that one and only one such candidate (IDM) is needed and its properties are constrained 
accordingly. This is of course not a fundamental principle and the possibility of multipartite 
dark matter should not be ignored. In this Letter we study its impact on the conventional 
picture of IDM physics, such as that of supersymmetry, using a simple generic scenario of 
two dark-matter candidates (2DM), one a fermion singlet (neutralino) and the other a scalar 
singlet. Our conclusions are broadly applicable to any 2DM model. 

Model The simplest way to have at least two DM candidates is to append the SM 
with the exactly conserved discrete symmetry Z2 x As pointed out in Ref. this 
may be realized naturally in the framework of = 2 supersymmetry. Alternatively, if 
the SM is extended to include an exactly conserved Z2 symmetry without supersymmetry, 
then the supersymmetric version of this extension will have Z2 x Z2, as in Refs. [Sj and [j]. 
To explore generically the impact of such a scenario, we first observe that the details of the 
specific model are mostly irrelevant, as far as the relic abundance, and the direct and indirect 
detection of dark matter are concerned, except for the masses of the two DM candidates and 
their interactions with the SM particles and with each other. This is because the relevant 
processes are either elastic scattering at almost zero momentum transfer or annihilation at 
rest. 

Specifically we add two new fields which are singlets under the SM gauge group: a new 
fermion x and a new scalar S. Under Z2 x Z2, x ~ (~) +) a-^id S ~ (+, — ), whereas all SM 
particles are (+,+). This means that {S) = is required. In a complete theory such as 
that of Ref. [3|, there may also be (— , — ) particles. For simplicity we assume that all such 
particles are heavy enough to decay into x ^"^^ <S- If not, we would then have to consider 
three coexisting DM candidates. 



3 



The Lagrangian of our generic 2DM model is given by 

1^ = ^SM + + ^DM + ^int, (l) 

where Csm denotes the usual SM Lagrangian, and 



1, ....... A3,,trr- , ^4 



A„i = -X^H^HSS + -j-H^Hxx + i^XXSS, (2) 

where H is the SM Higgs doublet. After the electroweak symmetry is spontaneously broken, 
H = {v + K)j\pl with V = 246 GeV and the masses of x ^^id S are given by = 
mi — X^v'^/2A and m| = + X^v^ 12. The various effective interaction terms, relevant to 
our discussion, are 

^.xx = 9MX. C,,,. = %hhxx, 

1 1 
^hss = ^9sv hSS, Chhss = -^9s hhSS, 

C^xss = ^XXSS, (3) 

where we have introduced the dimensionless couplings = X^v/A, gs = X2, and g^s = 
X4V/2A. Note that this is not meant to be an effective theory in powers of 1/A for all 
processes. It is applicable only to DM-nucleus elastic scattering (with almost zero momentum 
transfer) and DM annihilation at rest. 

As an example of how the effective couplings of Eq. (3) may be generated in a complete 
model, let us consider Ref. where a second pair of scalar superfields (77°, r/f ) and (r/^, 77°) 
are added, which are odd under a new Z2, whereas the usual and (02)02) 

MSSM are even. Together with the conventional R parity, we then have an exactly conserved 
Z2 X Z2 symmetry. Consider now the interaction 

IgyBi^iV'i - ^%°) + IgyBiM - M). (4) 

where B is the U{1)y gaugino. We may thus identify the x of our generic model with B, and 
S{h) with a linear combination of the real parts of 7?i2(0i2)- The effective XX'S'S and hhxx 
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interactions are then generated from the exchange of the rj and higgsinos, respectively. 
Assuming the masses of these higgisinos to be comparable to and m^, the effective 
couplings and g^s are not necessarily very much suppressed. This allows us to consider 
three characteristic scenarios, as depicted in Fig. [H For definiteness, we consider > ms 
in this analysis, but our conclusions are mostly the same if we switch them around. 

XX 
\ 

\ / 

SM 

(c) 

FIG. 1: Possible annihiliation scenarios in the 2DM model where the (red) arrow line denotes the 
DM annihilation. 

Scenario A [Fig. 1(a)]: g^ids 7^ but g^s = 0; both x and S can annihilate into SM 
particles but they do not interact with each other. Scenario B [Fig. 1(b)]: g^ = but 
95,9x3 7^ 0; X can only annihilate into S, after which S will annihilate into SM particles. 
All of x's information is hidden behind S, hence it will be superdark and cannot be observed 
directly. It may be revealed nevertheless if apparent discrepancies occur among different 
experiments. Scenario C [Fig. 1(c)]: gx,9s,9xs 7^ 0; x can annihilate into both S and SM 
particles, after which S wih annihilate into SM particles. The special case of g^s ^ 9x^9s 
is of particular interest. Here XX will annihilate predominantly into iSiS, resulting in a much 
smaller x I'elic abundance, thereby relaxing the constraints on its parameter space, which 
may be identified with that of the MSSM. Scenarios A and B are of course just two special 
limits of C, but they have qualitatively different predictions on the direct and indirect search 
experiments of dark matter, as shown below. 

Observational Constraints If two DM candidates coexist, the usual observational 
constraints also apply, but with modifications. 

(i) Relic abundance: Since both DM candidates contribute to the relic abundance, they 
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must add up to account for the current observation: 

n^h'^ + VLsh'^ = VLcDMh'^ = 0.110 ± 0.013. (5) 

It is well-known that the relic density of each DM species is approximately given by Qih"^ ^ 
(0.1 pb)/ {(jv)-, where (cv)- is the thermally averaged product of its annihilation cross section 
with its velocity. Using Eq. ([5]), we then obtain 



{av)^ ~ pb. (6) 



(a) Halo density profile: For simplicity, we assume the two DM candidates to have the 
same density profile and use that given by Navarro, Frenck and White (NFW) in our 
analysis. (It is of course straightforward to extend our results to other density profiles.) In 
the 2DM model, the dark-matter mass density profile of the galactic halo is thus given by 

p(r) = '-2^ , + '-^ (7) 

(r/rj (1 -I- r/rj (r/rj (1 + r/rj 

where Vc = 20.0 kpc and po is adjusted to reproduce the local halo density at the Earth 
position. Here, represents the fraction of the mass density of the ith dark matter in our 
local dark-matter halo as well as in the Universe, i.e. 

where pi is the local density of the ith DM and = 1. For our 2DM model, we obtain 

(Hi) Direct search: Assuming that DM is the dominant component of the halo of our 
galaxy, it is expected that a certain number of these weakly interacting massive particles 
(WIMPs) will cross the Earth at a reasonable rate and be detected by measuring the energy 
deposited in a low-background detector through the scattering of a WIMP with a nucleus 
of the detector. So far most experimental limits of this direct detection are given in terms 
of the cross section per nucleon under the IDM hypothesis. The event rate per unit time 
per nucleon is given by 

i?-5^n.(a), = 5^^(a),, (10) 
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where rij is the local number density of the ith DM and (a)- is the ith DM-nucleon elastic 
scattering cross section which is averaged over the relative DM velocity with respect to 
the detector. The measured experimental rate in the IDM case is given by -Rexp ~ PoC"o/'^o 
where cxo denotes the "zero-momentum-transfer" cross section of DM-nucleon scattering and 
mo is the DM mass. The current direct-search limit implies R < Rexp, i-e. 

— (^x^f ^SAf < — , (11) 

ms mo 

where c^AficrsM) denotes the scattering cross section of with a nucleon Af. Although 
the experimental sensitivities and limits are often described in terms of the dark-matter 
elastic scattering with a single nucleon, one should keep in mind that nuclear form factors 
may need to be taken into account. 

In Scenario B, there is no scattering of x with the nucleon, hence the limit in Eq. ffTTj) 
becomes ctsjv < <^o/^s, i-e- bounds from direct detection become weaker in this case. On 
the other hand, if dark matter is observed in direct-detection experiments, the DM-neucleon 
cross section may be understimated by a factor of l/eg- 

(iv) Indirect gamma-ray search: The relic dark matter may collect and become gravi- 
tationally bound to the center of the galaxy, the center of the Sun and the center of the 
Earth. If this happens, then a variety of indirect dark-matter detection opportunities arise. 
In particular, the measurement of secondary particles coming from dark-matter annihilation 
in the halo of the galaxy will help to decipher the nature of dark matter. Efforts to detect 
the annihilation products of dark-matter particles in the form of gamma rays, antimatter 
and neutrinos are collectively known as indirect detection. Of these, the observation through 
gamma rays is the simplest and most robust. The diffusion gamma-ray spectrum is given 
by 

d^ 2 d^x 2 d^s 



_ 2^^ , 2^^ flO^ 

dE^ ''dE^ ^dE^' ^ ' 

where d^i/dE^ {i = x, <S) is the differential gamma-ray flux along a direction that forms an 
angle ip with respect to the direction of the galactic center: 

2 

dl. (13) 



d^i dN^ , . 

{av) 



dE^ dE^' " Aiim^ [ (r/r,) (1 + r/r,)' . 



Po 
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The integral is performed along the line of sight. All annihilation channels of the ith DM 
are summed, and dN^/dE^ is the differential gamma spectrum per annihilation coming from 
the decay of annihilation products. 

Consider the special case of = mg = ruQ. After some simple algebra, one can show 
that 



d^ dK 



dEy dE^ ^ 47rmg 



Po 



-\ 2 



dl, (14) 



(r/rj (1 + r/rcf_ 

where we have used the fact that dN^/dE^ is almost the same for most of the final states. 
The integrated flux of the 2DM model is of the same order as that of the IDM model. 

(v) Collider search: Since Qh"^ oc 1/ (erf), the requirement of the correct relic density 
{QcDuh"^ ~ 0.1) implies that DM annihilation was efficient in the early Universe. It also 
suggests efficient annihilation now, implying large indirect detection rates, as well as efficient 
scattering now, implying large direct detection rates. The sum rule, cf. Eq. (jS]), means that 
the DM annihilation of each individual candidate has to be more efficient than that of 
the IDM case. Hence larger cross sections of DM production are expected at the collider. 
The smaller the fraction e^, the easier is the detection. In our simplistic case where the 
two DM candidates interact only with the SM Higgs boson, the vector-boson-fusion process 
qq qqVV qqh, with the subsequent decay h — > xx/SS, provides the most promising 
collider signature of the model ^] when ruh > m^/s- 

2DM Implications We first study the cosmological implications of either x or 5 as 
the sole source of dark matter. In Fig. [2|^a) and (b) we present the correlations between 
the effective coupling and the DM mass j?! , which is derived from WMAP data. The black- 
solid (red-dashed, blue-dotted) curve denotes flih"^ ~ 0.1 (0.05, 0.01), respectively. In the 
region below the black-solid curve the dark matter is overproduced. Fig. [2|c) shows the 
spin-independent cross section of DM-nucleon scattering for x (black-solid) and S (black- 
dashed). Current CDMS limit and projected sensitivity of CDMS2007 are also plotted. 
Using Eq. ffTTj) . we then derive a realistic bound on the 2DM model . For a large range of the 
DM mass, o"o/mj is almost a constant, e.g. (Jo/rrii ~ 2 x lO^^pb/GeV (2 x 10^^°pb/GeV) 
for the current CDMS data (projected CDMS2007 sensitivity). In Fig. EJ^d) we present 
the allowed parameter space of Scenario A in the 2DM model in the {ms, m^) plane for 
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FIG. 2: (a) and (b) show the correlations between the couphng and the mass of a single DM 
candidate as determined by the WMAP data: (a) for x (b) for S. (c) shows the spin- 
independent cross sections of DM-nucleon scattering in the IDM model together with the CDMS 
limit and future projected sensitivities of CDMS2007. We choose ruh = 200 GeV throughout in 
this work, (d) shows the allowed {m^, mg) parameter space of Scenario A in the 2DM model. 



= es = 0.5. In Scenario B, the limits only depend on S and the bounds become weaker. 

A promising way to tell if there are two coexisting DM candidates, assuming that they 
are very different in mass, is through indirect gamma-ray observations. The overlap of 
the two distributions might change the line shape of the gamma-ray distribution which 
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FIG. 3: Predicted gamma-ray spectra in the 2DM model for e-^ = es = 0.5. The predicted gamma 
flux is from a AO = 10~^ srad region around the direction of the galactic center, assuming the 
NFW halo profile (with a boost factor as indicated in the figure). For comparison we also show the 
scaled gamma-ray distribution in the IDM case. EGRET and HESS observations are also shown 
here for comparison. 

is distinguishable from that of the IDM case. In Fig. [3] this is illustated by showing the 
predicted fluxes from a /S.Vt = 10^^ srad region around the direction of the galactic center 
together with the existing EGRET [9] and HESS 10| observations in the same sky direction. 
We adopt the NFW density profile for the DM in our galaxy (J x AQ ~ 1 for AQ = 
10~^ srad) and allow the flux to be scaled by a "boost factor". For demonstration we choose 
= es = 0.5, = 260 GeV, and = 60 GeV. Clearly, the resulting gamma-ray flux 
distribution from the overlap of 2DM distribution is significantly different from that of the 
IDM model, which can be probed by the GLAST experiment HI]. The gamma-ray spectra 
can also be used to distinguish Scenario A from B of the 2DM model. Unfortunately, it is 
difficult to observe a shape change if — ms is small. On the other hand, they may be 
discriminated at the Large Hadron Collider (LHC) because in Scenario A, both x and S are 
produced; whereas in Scenario B, only S is. A discrepancy between relic abundance and 
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FIG. 4: Relation between g^. and when the XX ~^ 'SS annihilation mode is open. For illustra- 
tion, we choose e-^{es) = 0.9(0.1) and ms = 120 GeV. 

LHC production may reveal Scenario B. 

Consider the special case {g^s ^ fi'x' 9s) Scenario C, where the new annihilation 
channel xX ~^ <S<S opens. This case is very interesting because it has a crucial impact 
on the conventional supersymmetric DM model. For example, the lightest neutralino is 
a well-motivated dark-matter candidate, but its relic abundance is typically too large, or 
equivalently, its annihilation rate is too small. The WMAP data thus impose very tight 
constraints on the parameter space of the MSSM. But those constraints can be relaxed if 
there exists an additional DM candidate which opens up a new annihilation channel for 
the neutralino. For illustration, we choose g^s = (0.5, 0.55, 0.6) and nis = 120 GeV with 
e^{es) = 0.9(0.1) in the 2DM model. Using Fig. Mh), we then fix gg = 0.114. The (black) 
dotted curve in Fig. H] denotes Q^h"^ = 0.1 in the IDM model and the region below it will 
exceed the relic abundance. After including the new annihilaton channel xX ~^ more of 
the parameter space is reclaimed. Increasing g^s will open up even more parameter space. 

Conclusion In this Letter we presented a simple generic model of two coexisting dark- 
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matter candidates. We discussed its three characteristic annihilation scenarios and its impact 
on the observational constraints of dark matter. We note that the cosmic gamma-ray ob- 
servation is a good probe for confirming the 2DM model. We also demonstrate that with 
a second dark-matter candidate, the usual severe constraints on the parameter space of the 
MSSM can be relaxed. More detailed studies of this new idea of multipartite dark matter 
are forthcoming. 
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